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Abstract 



To any compact Riemann surface of genus g one may assign a principally polarized abelian 
variety of dimension g, the Jacobian of the Riemann surface. The Jacobian is a complex 

■ torus and we call a Gram matrix of the lattice of a Jacobian a period Gram matrix. We 
, give upper and lower bounds for all entries of the period Gram matrix with respect to a 

suitable homology basis. These bounds depend on the geometry of the cut locus of non- 
separating simple closed geodesies (scg). We first present a theoretical approach that relies 
on the premise that these cut loci can be calculated. Then we give an example where our 
, , upper bound is sharp. Finally we give practical estimates based on the geometry of Q-pieces 

■ that contain a canonical homology basis. The methods developed here have been applied to 
\ surfaces that contain small non-separating scgs in [ BMMSj . 

O ; Mathematics Subject Classifications: 14H40, 14H42, 30F15 and 30F45 

CN ■ 1 Introduction 

Let S be a hyperbolic Riemann surface of genus g > 2 and let A = (ai)j^]^ 2g ^ canonical 
K / , homology basis of closed cycles on the surface. This basis is assumed to be given in the following 

; way: 

.5^. A = (a;l,a.^(l),•••,"g,ar(g))■ 

Here each Oi is a simple closed curve and the curves are paired, such that for each there exists 
exactly one aT-(j) G ictj)j^i 2g ^^^^ intersects ai in exactly one point and there are no other 
intersection points. In the vector space of real harmonic forms on S let (o"fc)fc=i 2g ^ '^^'^^ 
basis for {ai)i^i 2g defined by 



The 'period Gram matrix Qs of a compact Riemann surface (R.S.) S of genus g is the Gram 
matrix 
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QS = i{(^i,(^j))ij = 1...2g = 

\S / i,j=1...2g 

This period matrix Qs defines a complex torus, the Jacobian or Jacobian variety J{S) of the 
Riemaim surface S (see [FK| . chapter III). We denote by 

s 

the energy of at (over S). As Qs is a Gram matrix E{ai) is also the squared norm of a vector 
Vi in the lattice of the Jacobian. By Riemann's period relations the Jacobian is a principally 
polarized abelian variety (PPAV) (see |BL| . section 4.1 for a definition). The lattices of PPAVs 
are exactly the symplectic lattices ( |BS| ). The Schottky problem is to characterize the Jacobians 
among the PPAVs. 

Buser and Sarnak ( |BS| ) approached this problem by means of a geometric invariant, the squared 
norm of the shortest non-zero vector, in the lattice of a Jacobian of a Riemann surface of genus 
g is bounded from above by log(45f), whereas it can be of order g for the lattice of a PPAV of 
dimension g. More insight has been obtained recently into the connection between the global 
geometry of a compact Riemann surface of genus g and the geometry of its Jacobian. The 
log(5()-bound on the squared norm of a lattice vector of the Jacobian has been further extended 
in [BPS) to almost g linear independent vectors. In |Mu3| it was shown that in the case of a 
hyperelliptic surface the squared length of the shortest lattice vector is bounded from above by 
a constant independent of the genus. 

In this article we examine the connection between the concrete geometry of a compact Riemann 
surface and the geometry of its Jacobian. This approach has been taken in special cases, for 
example when the Riemann surface is a real algebraic curve. Here algorithms exist to calculate 
the period matrix (see ^BSij . ^Hij or [Se]). but these methods are limited to these cases. Though 
the estimates obtained in this article are not sharp in general, the methods apply to any Riemann 
surface. 

Here we give upper and lower bounds for all entries of the period Gram matrix given with respect 
to a canonical homology basis. To this end we have to estimate the energy of the dual harmonic 
forms on the surface. We find upper bounds for the energy of these forms by estimating the 
capacity of hyperbolic tubes. The boundary of such a tube is provided by the cut locus of a 
cycle in the homology basis or a related scg. We then use projections of vector fields onto curves 
to obtain a lower bound. This theoretical approach relies on the premise that these cut loci 
can be calculated. We furthermore provide two examples for our approach. One that shows the 
limitations of the method and one where the upper bound is sharp. Finally we present practi- 
cal estimates based on the geometry of surfaces of signature (1,1), or Q-pieces, containing the 
paired curves of the homology basis. Under this condition the cut loci can be - at least partially 
- calculated. The methods developed here have been applied to surfaces that contain small scgs 
in |BMMS| . In this case more refined estimates are obtained. 

This article is structured in the following way. We will first introduce the necessary tools and 
definitions in section 2. Then we present the theoretical approach in section 3. Finally we provide 
practical upper and lower bounds based on the geometry of Q-pieces in section 4. 
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2 Preliminaries 



Many calculations presented in the following sections rely on the embedding of topological tubes 
around scgs of Riemann surfaces into hyperbolic cylinders and subsequent approximations and 
calculations in Fermi coordinates. We will present all these concepts in this section. 

The Poincare model of the hyperbolic plane is the following subset of the complex plane C: 

M = {z = x + iyeC\y>0} 

with the hyperbolic metric 



y 

Fermi coordinates i/j^ with baseline 77 and base point p are defined in the following way. The 
Fermi coordinates are a bijective parametrization of H 

V',, : ^ H, V'r; : {t, s) tprjit, s), 

where ijjr](0,0) = p. Each point q = tprjitjs) G H can be reached in the following way. Starting 
from the base point p we first move along r/ the directed distance t to ipr]{i, 0)- There is a unique 
geodesic, v, intersecting rj perpendicularly in ^^(t,0). From tp^{t,0) we now move along u the 
directed distance s to ^rj{t,s). 
The parametrization : — >■ H 

cosh(s) 

are the Fermi coordinates with baseline {iy \ y G ^+} and base point i. Fermi coordinates for 
any geodesic in EI can be obtained by conjugation of s) with the suitable isometries. 
A hyperbolic cylinder C or shortly cylinder is a set that is isometric to the set 

{i^{t, s) I {t, s) G [0, a] X [61, 62]} mod {^^(0, s) = V'(a, s) | s G [61, 62]}, 

with the induced metric from H. The baseline of C is the shortest scg in C of length a. 
Consider a cylinder C. Let A C C be a set and F G Lip(^) a Lipschitz function on the closure 
of A. Let G be the metric tensor with respect to the Fermi coordinates. Then the energy of F 
on A, Ea{F) is given by 



Ea{.F)= J J \\D{Fo^)\\l_, v/d<G). 



V>-i(^) 

Using Fermi coordinates, we obtain for the energy of F on A, Ea{F), with F otp = f: 

The capacity of an annulus A G C, cap (A) is given by 

cap(^) = mi{EA{F) \ {F G Lip(^) | F\d^^A) = O.i^la^CA) = 1}}- 
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Let S be a compact R.S. and {cti)^^-^ 2g ^ canonical liomology basis on S. The collar C{r]) of a 
simple closed geodesic rj, is defined by 

C(r/) = {x £ S \ dist(x,r/) < w} . 

Here w is the supremum of all u, such that the geodesic arcs of length uj emanating perpendic- 
ularly from r] are pairwise disjoint. The closure C{t]) of a collar is a cylinder of constant width. 
For a given Oj, a^(^i-j is the unique scg in the canonical homology basis that intersects Oj. In 
|Mu3| it was shown that the capacity of any cylinder C with baseline a^(j) contained in C(a^(j)) 
provides an upper bound for the energy E{(Ti) of the dual harmonic form cjj. If u)T-(j) is the width 
of C(aT-(j)), then 

E{ai) = [aiA *ai < ^^^^ ^ = cap(C(a,(,))) < capiC). (2) 

S ^-2arcsin(^^;^j^(^j 

It is noteworthy that we can deduce this inequality by finding the minimal function for the 
capacity problem. We obtain this function by solving the Euler-Lagrange equation (see |Ge) p. 
152-154) applied to equation ([T]). For a more detailed account see ^Mulj . p. 23-27. 



3 Theoretical estimates for the period Gram matrix 

Let 5" be a Riemann surface of genus g and A a canonical homology basis. Let Qs be the period 
Gram matrix of S with respect to A, 

..=(/.-.) 

\'S' / i,j=1...2g 

From equation ([2]) we obtain an upper bound on the diagonal entries of Qs by evaluating the 
capacity of a cylinder around the elements of the canonical homology basis. This approach can 
be expanded to obtain better estimates of these entries and finally of all entries of Qs- The 
approach relies on the premise that the cut locus of a given simple closed geodesic on a Riemann 
surface can be - at least partially - calculated. 



3.1 Estimates for the diagonal entries oi Qs 

We first show how to obtain upper and lower bounds on the diagonal entries of Qs- Consider 
WLOG E{ai) = qii- Let T{a2) C S he a topological tube which is obtained by a continuous 
deformation of the collar C(a2) in S- Using the same arguments that are used in iBSj, p. 36 for 
hyperbolic cylinders, we see that the capacity of such a tube T{a2) gives an upper bound for the 
energy of ui, E{ai): 

E{ai) < cap(r(Q2)). 

We will use the tube that we obtain by cutting open S along the cut locus CL{a2) of a2- The cut 
locus of a subset X d S, CL[X) is defined in the following way, where ^a,b denotes a geodesic 
arc connecting the points a and h: 

CL{X) ■= {y e S \ ^'^x,y,lx',y-,lx,y / lx',y-, with X, x' € X and dist(x,y) = l{-ix,y) = l{lx',y)}- 
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For more information about the cut locus, see |Ba| . We denote by Sx the surface, which we 
obtain by cutting open S along CL[X). For a set X d S, set 

Zr{X) = {xeS \ dist(x,X) < r}. 

Let U he a. set of disjoint simple closed geodesies (7i)i=i,...,n then for sufficiently small r, Zr{U) 
consists of disjoint cylinders around these geodesies. We obtain CL[U) by letting r grow con- 
tinuously until Zr{U) self-intersects. We stop the expansion in the points of intersection, but 
continue expanding the rest of the set, until the process halts. The points of intersection then 
form CL{U). It follows from this process that the surface Su, which we obtain by cutting open 
S along CL{U) can be retracted onto the union of small cylinders around the (7j)j=i,...,n- If 
f7 = 7, then Su can be embedded into an sufficiently large cylinder C around 7. 
Consider an embedding of 8^2 = 5*2 in a cylinder C, which, by abuse of notation, we also call 
82- This is shown in Figure{l\ The boundary 882 of 82 d C consists of the two connected com- 
ponents 8182 and 8282, which are piecewise geodesic (see |Ba| ). Fixing a base point x G 82 C C, 




Figure 1: Embedding of ^2 = 8a2 in a cylinder around 02 • 

we can construct a primitive Fi of ai by integrating ai along paths starting from the base point 
X. As J 0"! = 0, the value of the integral is independent of the chosen path in 5*2. Therefore 

a2 

there exists a primitive Fi of ai on 5*2 C C. Furthermore Fi is a real harmonic function, as o"i 
is a real harmonic 1-form. We remind here that the value of the integral of a ai over a closed 
curve depends only on the homology class of the curve. Especially the value of the integral is 
the same for two curves in the same free homotopy class. 

The conditions on the homology basis A imply the following boundary conditions for Fi. For 
each point pi on the boundary 882 C C, there exists a point p2 that maps to the same point p 
on 8 as pi, such that 

Fi{p2) - Fiipi) = or Fi{p2) - Fiipi) = 1. 
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We color pi and P2 blue in the first case and red in the second case and call such a decomposition 
a red-blue decomposition of the cut locus. Such a decomposition is shown in Figure{l\ 
We denote by CL^'"^(a2) and CL''''^°'(a2) the blue and red parts of CL{a2), both in S and 82- 
We have that 

CL{a2) = CL^^''%a2) U CU^'^{a2). 

For the red-blue decomposition that is obtained via the cut locus CL(a2), the following holds. 
If pi and p2 are blue, then pi and p2 lie on the same side of the boundary dS2- If pi and p2 are 
red then they lie on different sides of dS2- This follows from the relationship of the canonical 
1-forms with the intersection number of curves (see |FK| . chapter III). There are points on 882 
that are both red and blue. However the number of those points is finite. 

We now connect the endpoints of two corresponding opposite red boundary segments in the red- 
blue composition of ^2 C C with differentiable curves, such that these curves do not mutually 
intersect. Then the curves, together with the boundary segments of S2, enclose a subset of 82- 
We denote the union of all enclosed areas that can be obtained in this way by 5f ^. We set: 

= 52\55'^ therefore S2 = S^"^ U 5^'"^ 

We also call such a decomposition of S2 a red-blue decompositon. Note that the connected red 
and blue parts on each boundary of CL{a2) alternate and that the boundaries of S^'"*^ contain 
the blue parts of CL{a2)- By construction of the red and blue parts of 5*2 do also alternate. 
The connected subsets of CL*'"^(a2) are paired via the red-blue decomposition. This induces a 
pairing of the connected areas of Two such pairs are isometric. 



Upper bound 



The capacity of 5*2 C C, cap(S'2) provides an upper bound on the energy of ai. We obtain 
an upper bound on cap(52) by evaluating the energy of any test function Fn, which is a Lip- 
schitz function on 5*2 and satisfies the boundary conditions of the capacity problem (see |GT| ): 

E{Fu) > cs.v{S2) > E{ai) = qu. 

In |Mu2| Theorem 4.2. we obtain general upper and lower bounds on the capacity of annuli 
^ on a cylinder of constant curvature. Here an annulus is a region that can be obtained by a 
continuous deformation of the cylinder itself. The idea for the upper bound is to construct a test 
function in the following way. We adapt the harmonic function that solves the capacity problem 
for cylinders of constant width (see ([2])) to the boundary using the parametrization of a cylinder 
in Fermi coordinates. If the annulus ^ C C is given in Fermi coordinates by 

A = V{(t,s) I [ai(t),a2(t)],t G [0,/(q2)]}, 
where ai(-) and a2(-) are piecewise differentiable functions with respect to t, then it follows: 

Theorem 3.1 There exists a Lipschitz function F G Lip(A), such that for H{s) = 2 arctan(exp(s)) 
and qi{t) = ^^^^\so=a,{t) ■ a'i{t) /or i G {1, 2}, we have: 

— — 1- — --— dt = EiF) > cap(yl) > / -— — — - — --- dt 
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As 5*2 C C is an annulus that satisfies the conditions of Theorem 13.11 we can apply this the- 
orem to obtain an estimate of the capacity cap(S'2). It is noteworthy that the lower bound is 
obtained by finding the optimal function that minimizes the last integral in equation ([1]). This 
integral is equal to the energy of the projection of the gradient of a function. Here the gradient 
is projected onto the geodesies intersecting the baseline of C perpendicularly. 
We now shortly describe a further modification that should lead to an improved upper bound. 
Though we have not applied this modification in our practical estimates, we think that it leads 
to a test function Fu that is close to the harmonic function Fi . This modification can be applied 
if a connected subset of 82^^ is large. On the boundary diS^ , where our test function Fn 
equals 1, we move along the boundary line and decrease the boundary values until they reach 
^. After a plateau, where the boundary values are constant we increase the values in the same 
fashion when towards the other end of diS^ . On the opposite boundary 828^ , we do the opposite. 
First we increase the values to ^. After a plateau we decrease the boundary values again to 0. 
There is a paired subset S'^ of 82^^ isometric to . Transplanting Fn from S'^ via the isometry 
onto S'^, we obtain a function Fu that satisfies the boundary conditions induced by the red-blue 
decompositon of Fi. However, in this case we have to assure that the 1-form dF\t extends to 
a closed 1-form on S. If this condition is satisfied, we have E{Fit) > E{ai) due to the energy 
minimizing property of ai. The energy E{Fit) of our new test function should provide a better 
upper bound for E{ai) than cap(S'2). 

Lower hound 

We obtain a lower bound on qn = E{Fi) in the following way. Consider the set 82^'^. Re- 
member that in each connected subset of 82^'^ there are boundary points pi and p2 on opposite 
sides, such that Fi{p2) — Fi{pi) = 1. We have: 

E{ai) = E{Fi)> j WDFiWl 

Qred 

Let / be a disjoint union of intervals in M and 

cp: I x[ai,a2]^ 81""^, if : {t, s) ^ ip{t, s) 
a bijective function that parametrizes 82'^'^ in the following way: 

(^(/ X {ai}) = 8^"'^ n 91^2 and ip{I x {02}) = 8^"'^ n 8282 
and for a fixed I, ^{{c} x [01,02]) is a differentiable curve in 82'^'^, such that 

Fi{ip{c,a2)) - Fi (93(0, 01)) = 1. 
Denote by F'l the set of functions 

Ti = {f:8'2'"'^m\feUpi8'2"') and /((^(c, 02)) - /Mc, oi)) = 1 V c G /}. 
We obtain a lower bound on qn = E{ai) = E[Fi) if we can find a function /i, such that 

j ||i^/i||i= mm j WDfWl (3) 

Qred Qred 
"^2 '^2 
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We call this problem the free boundary problem for 82'^'^. 

Though this problem is quite interesting in its own right, we could not find an explicit solution. 
To obtain an explicit result, we construct another lower bound based on projection of tangent 
vectors onto curves. For (p{c, a) e 5f denote by : T^{S^^'^) {X ■ | A e M} the 

orthogonal projection of a tangent vector in x onto the subspace spanned by ^'^g^'"^ ■ Then we 
have: 

E{F,)> I \\DF,\\l> I \\p^{DF^)\\l>n^n I \\p^{Df)\\l= J \\p^{Dh)\\l (4) 

Qred Qred Qred Qred 

J2 '-'2 '-'2 '-'2 

Here fi is a function that realizes the minimum. We have / HDFiHl = J \\py:,{DFi)\\2, if 

Sred Qred 
2 '^2 

and only if in every point ip{t, s) = x G ip{t, •) is orthogonal to the level set of Fi passing 
through X. Note that the problem of finding the function fi is in general easier than finding the 
function Fi or fi. We will apply these ideas to Q-pieces in the following section. To this end we 
will use results from the calculus of variations. 
In total we have the following estimates for qn = E{ai): 

cap{S^,) > E{ai) = E{F,) > min [ > min [ \\p^iDf)\\l 

Qred Qred 
■^2 "-^2 

If the geometry of the cut locus CL(a2) is known, we obtain an explicit estimate for cap(5Q,2) 
from Theorem 13.11 

Note that the test function whose energy provides our upper bound has positive energy on S'g^"*^, 
whereas the energy is zero on 82^'^ in the estimate providing the lower bound. This is one 
reasons why the upper bound differs from the lower bound. Another difference is due to the use 
of the projection along lines in the construction of the lower bound. In section 4, we will apply 
these methods to a decomposition of the Riemann surface, where the elements of the canonical 
homology basis are contained in Q-pieces. Here we will see these two effects explicitely. 

3.2 Estimates for the non-diagonal entries of Qs 

We now show, how we can estimate the remaining entries of the period Gram matrix Qs- For 
i ^ j, we have, as J • A *• is a scalar product: 
S 

h,\ < ^{E{ai) + E{a,)), (5) 
Qij = ^{E{ai + aj) - E{ai) - E{aj)) and (6) 
qij = ^{E{ai) + Eiaj)-E{a,-aj)). (7) 

We have shown how to find upper and lower bounds on E{ai) and E(aj). We obtain a direct 
estimate on qij from equation ([5]). However to obtain a sharp estimate, both E{(Ji) and E{aj) 
must be small. We will show how to obtain better estimates for qij from the following two 
equations. Here we obtain an estimate for qij, if we can find upper and lower bounds on either 
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E{ai + CTj) or E{ai - aj). 

Consider the harmonic 1-forms wi and ^2 that satisfy the following equations on the cycles: 

J uji = Sik + 5jk for all ke{l,..,2g}. (8) 

y ^2 = Sik-Sjk for all ke{l,..,2g}. (9) 

It follows easily from the uniqueness of the harmonic 1-form that satisfies these integral conditions 
that (Tj + cjj is the unique 1-form that minimizes the energy among all closed 1-forms that satisfy 
equation ([8]) and that — aj is the unique 1-form that minimizes the energy among all closed 
1-forms that satisfy equation Q. 

There is a geodesic a in the homotopy class of either 0^(4) • Q;^(j) or aT-(j) (aT-(j) )~"'^ which is a simple 
closed curve. On Sa the primitive of either ai+aj or —aj induces a red-blue decomposition of 
Sa (see for example (T2 + (T4 on Sai-a^ in Figure{l^. Therefore we can proceed as in the previous 
subsection to obtain estimates on CTj + aj or fjj — aj on Sa and hence for qij using equation ([6]) 
and ©. 

As it is in general difficult to find a, we will present this approach only in the case, where 
aj = We present these estimates in Case 1. If aj 7^ ar(i)) we will present an alternative 

approach in Case 2. We will make use of both methods in section 4.2. 

Case 1: Estimates for Qi^^^i^ 

Consider WLOG qi2- Consider the simple closed geodesic Q12 in the free homotopy class of 
ai2 = Oil 02"^- Due to the relationships in equation ([8]), ai+a2 has a primitive F12 on Sai2 = S12 
that induces a red-blue decompositon of 812- We embed S12 into a cylinder C. We also denote 
the embedded surface by 512- 

As in the previous subsection, the capacity cap(S'i2) provides an upper bound for E{ai + o"2): 

cap(S'i2) > E{ai + (T2). 

We obtain a lower bound for E{Fi2) = E{ai + (T2) from the red-blue decomposition induced by 
F12 on on 512. Here we can use the same methods as for E{ai) on 5*2 presented in the previous 
subsection. 

With the estimates for E{ai +1T2), E{ai) and E{(T2) applied to equation ([6]), we obtain an upper 
and lower bound on qi2- 

Case 2: Estimates for qij,j 7^ r(i) 

In this case ai and aj do not intersect. Consider WLOG 513. For aT-(i) = 02 and 0^(3) = 04 
consider ua4 = 5'24- 5*24 consists out of two connected parts. Let S^^ C be the part that 
contains 02 and let '^"4 be the part that contains a^. We embed into a cylinder Ci 

around 02 and 524 into a cylinder C3 around and denote the embedded surfaces by the same 
name. Due to the relationships in equation ([8]), ui + (T3 has a primitive on both 524 C Ci and 
5'24 C C3. Such a decomposition is shown in Figure\^ 

For i G {1,3}, let Fi on S'24 be a function that satisfies boundary conditions for the capacity 
problem on S'24. Then these functions together define naturally a function F13 on S'24. By 
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Figure 2: Embedding of in a cylinder Cj around a^i^i) for i G {1,3}. 

smoothing F13 in an inner environment of the boundary of 5*24, we obtain a function /13 on 
S, whose derivative d/13 is a closed differential form that satisfies the same integral conditions 
on the cycles as cxi + CT3. Due to the energy- minimizing property of ai + (73 we have that 
E{fi3) > E{ai + (T3). Hence the sum of the capacities of S^^ and S'24 provides an upper bound 
for E{ai + (73): 

cap(5'24) + cap(S'|4) > E{ai + 0-3). 

We obtain a lower bound for E(ai + (73) by applying the methods used to obtain a lower bound 
on E{ai) on S2 presented in the previous subsection. Here we obtain estimates from the red-blue 
decompostions induced by the primitive Fi^ of cti + (T3 on S'24 ^1 ^^'^ ^24 ^3- '^^^ only 
difference is that we have some segments of the boundary, where the red-blue decomposition 
does not apply. Here we disregard these pieces in the construction of 5'24^'^ and §24^'^ • As these 
sets are disjoint we are looking for a certain function /13 that satisfies 

/isfe) - fisiPi) = Fn{p2) - Fn{pi) = 1 

for all points pi,p2 in 95'24^'^ or dS^^^'^- Let p^p be the projection of a vector field in the tangent 
space onto lines of a suitable parametrization 99 of S'24'^'^ and 524^"^. Let /13 is the function that 
minimizes the projected energy E{j)^{D-)). We then obtain with S'24 ~ S'24'^'^ U S'fl^'^: 

E{ai + a:i)>EsrJp^(Dh:,)). 

With the estimates for E{(Ti +C73), E{ai) and E{a3) applied to equation ([6]), we obtain an upper 
and lower bound on qi^: 

qi3 > ^ {Es^^{p^{Dfis)) - cap(SQ,2) - cap(S'Q,4)) < and 
gi3 < I (cMSL) + cap(5i4) - Esre,{p^{Dh)) - Esr.,{p^{Dh))) > 0. (10) 
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Here is the minimizing function corresponding to the primitve of on constructed 
analogously to fi in the previous subsection. Using this approach we can only obtain optimal 
estimates, if qi^ is close to zero. This is due to the fact that we do not have full information of 
the boundary values on our tubes S'24 and 6*24. This estimate is however better than the one 
obtained from equation ([5]). As was shown in |BMMS| , the value of qi^ is always close to zero, 
if Q2 and 04 are separated by a small scg 7. 

We now give two examples to demonstrate the weakness and strength of the method. We first 
show that the energy of a dual harmonic form can be lower than the capacity of a cylinder of 
even infinite length. The following example is due to Peter Buser. 




Figure 3: Building blocks for the surfaces N and L of genus g. 

Example 1: For comparison we briefly review here the example of the necklace surface given 
in [BSelJ. Let y be a Y-piece, a surface of signature (0,3). Let 7,77 and 77' be its boundary 
geodesies, such that r] and r]' have equal length. We paste two copies of Y along ij and r]' to 
obtain R of signature (1,2). The necklace surface of genus g, is obtained by pasting together 
g — 1 copies Ri, ...Rg-i of a building block R as shown in Figurel^ Here the free boundary of 
Rg-i is pasted along 71 of Ri to obtain a ring. By the collar lemma (see |BuJ. p. 106), each 7^ 
has a collar of width w-y, where 

> arcsinh ' 



^ \sinh(^) 

Let A = (ai, a.^^!), a^, aT-(g)) be a canonical homology basis, such that ai = 71 and «,-(!) = 02 
is a simple closed geodesic that intersects all 7^ exactly once. Let Qs be the corresponding period 
Gram matrix. We will examine the upper bound on the entry 522 = E{a2)- 

Following our method we have to embed 5*1 = Sa^ into a cylinder Ci and have to evaluate 
cap(5i). Now if an annulus A is contained in an annulus B then cap(A) > cap(i?). Therefore 
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cap(5'i) > cap(Ci). The capacity of the cyhnder Ci with basehne of length l{ai) = l{'y) of 
infinite width is not zero. This follows from equation ([2]). We have that 

cap(5i) > cap(Ci) = 

vr 

We now give another estimate for the energy of a2 with the help of a test form S2 ■ This approach 
applies only to this example. To this end consider the collar C(7j) of a 7j. On each C(7j) set 
S2 = dF2 , where dF2 is the real harmonic function that has value on one boundary of C^'ji) 

9-1 

and — j- on the other. We set S2 = on 5\ |J C(7i). Then S2 is arbitrarily close to a closed 

i=i 

form that satisfies the same conditions on the cycles of A as a2 and we have 

vr — 2 arcsin ( — t-, — ^ ) 

where is bounded from below by the collar lemma. Hence E{a2) is at most of order ^ and 
goes to zero, as g goes to infinity. Our upper bound, on the contrary, is always bigger than the 
constant This shows that there exist examples where our upper bound can not be of the 

right order. This might be due to the fact that the projection of CL^^^'^{ai) onto ai can attain 
almost the length of ai. Hence as 5"^'"^ is large, E{a2), the energy of a2 might be small. However 
it might be possible to obtain a better result by adjusting the boundary values on CL^'"'^(ai) as 
described in the discussion of the upper bound. 

Example 2: For our second example we construct a linear surface L of genus g. This ex- 
ample belongs to the class of M-curves described in (BSi| . In this construction we use Y-pieces 
Y, where the length of r] and rj' is large. We construct R from two copies of these Y-pieces as 
in the previous example. To construct a surface L of genus g, we paste together g — 2 copies 
i?2) ■■■Rg-i along the 7^. Then we paste two copies of Y , Yi and Yg on each side of R2 and Rg-i, 
respectively. We close these Y pieces by pasting them together along r] and r]' . We call the image 
oi rj in Yi r]i for i G {1,9} (see Figurel^. 

Let A = («!, a^(i) , ...,ag, ct^(^g)) be a canonical homology basis, such that ai = r]i and = 02 
the unique simple closed geodesic in Yi that intersects ai perpendicularly. Let be the cor- 
responding period matrix. We now show that in this case, the upper bound for 522 = E{a2) is 
optimal. Therefore we use the symmetries of the surface L. 

It is easy to see that CL(ai), the cut locus of ai in L consists of the scgs ul that intersect the 
geodesies 7^ perpendicularly and the scgs v\ C Ri that intersect the scgs r/j and 77^ perpendicularly 
(see Figure\^. We have: 

CL(air^ = K,...,i/^^} and CL{aif'^- = {ui...,vl^,}yjrjg. 

We now consider the following two isometrics (pi and 02 in Isom(L). 0i is the hyperelliptic 
involution that fixes the cut locus CL{ai) as a set, such that 4'i{rji) = r]'- for all i. (f)2 is the 
isometry that fixes CL(ai) as a set and all u'^ point-wise. Set (p = (pio (j)2- (p is the isometry that 
maps any point p in CL{aiY^'^ C to the corresponding point p' in the red-blue composition 
induced by a2- 

Consider now the primitive F2 of CJ2 on Sa^ = 5*1. We have that F2 o (p2 is a, harmonic function. 
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whose derivative d{F2 o (f)2) defines a 1-form on S. a'^ satisfies the same conditions on the 
cycles as 02- Due to the uniqueness of a"2 we have a'2 = <72 and we have -F2 o 02 = ^2- In the 
same way 1 — i<2 o is a harmonic form, whose derivative —d{F2 o (pi) defines a 1-form a!^ on 
5 that satisfies the same integral conditions on the cycles as CJ2. Hence we have a2 = o"2. By 
choosing an appropriate additive constant, we obtain: 

F20 (p2 = F2 and 1 — i<2 o 01 = F2. 

Now for any p on one side of CL{aiY'^'^ C 5i we have that F2{p) — F2{(l){p)) = 1. Together with 
the above equation we obtain: 

^2(0(p))=O and F2{p) = l. 

Hence the red parts of the boundary satisfy the conditions for the capacity problem. Now for 
large r/, we have from hyperbolic geometry that is arbitrarily small. Hence in this case we 

obtain 

q22 = E{a2) = cap(S'i) - e 
and our upper bound for a diagonal entry of Qs is sharp. 

4 Estimates for the period Gram matrix based on Q-pieces 

In this section we establish upper and lower bounds for all entries of the period Gram matrix 
based on the geometry of the Q-pieces. All formulas concerning the geometry of hyperbolic 
polygons in this section can be found in [BuJ, p. 454. 

Let 5 is a Riemann surface of genus g > 2 and A a canonical homology basis. Then the pairs 
((oj, a,-(j)))j=i^3^..23-i are contained in disjoint surfaces of signature (1, 1) or Q-pieces. If g = 2, 
then the Q-pieces intersect at the boundary, but this will not concern us here. Let (Qi)i=i,3,--,23-i 
be a set of Q-pieces, such that 

{ai,ar(i)) C Qi = Qr(i) for i G {1, 3, .., 2^- - 1}. 

We will present estimates on the entries of the period Gram matrix Qs using a partition of S 
that contains these pieces. To establish these estimates, the following information about the 
geometry of the Q-pieces (Qi)i=i,3,..,2g-i must be known. 

Let /3j be the boundary geodesic of Qi and let aiT-(j) C Qi be the simple closed geodesic in the 
free homotopy class of aj(aT-(j))~^. If we cut open Qi along one of these three geodesies, Qi 
decomposes into a Y-piece. Let tj g] — ^, ^] be the twist parameter (see |Bu) . p. 29) for the 
decomposition with Oj, where j £ {z, r(z), ir(i)}. To obtain practical estimates on all entries of 
the matrix Qs, we assume that the following parameters are known for each Qi: 

- the length of /3j 

- the length of Oj, a^i^i-^ and aj 

- the twist parameters tj for j G {i,T{i),iT{i)} 

A Q-piece Qi is already determined by a triplett {/3i,ak,tk), where k G {i, r(i), ir(i)} and these 
Fenchel-Nielsen coordinates are a common way to provide the data for a Q-piece. Therefore we 
will first show how to obtain the remaining information in the list above from such a description. 
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Then we will give a parametrization of trirectangles. The formulas developed here, will be used 
in the final subsection. All formulas in this section have been implemented in a MAPLE sheet, 
so that concrete estimates can be obtained from the above input. 

Let WLOG Qi = Qi- We assume that Qi is given in the Fenchel-Nielsen coordinates ai,ti). 







02' 





Figure 4: Two lifts of ai in the universal covering 
We assume furthermore that 

cosh(^) < cosh(^) + i. (11) 

Such a short ai exists indeed. This is stated in [Pa], Proposition 5.4. This assumption 
simplifies the calculation of the twist parameters t2 and ti2 ■ As we will see in the following this 
choice of a small ai also improves our estimates. 

In Qi there exists a unique shortest geodesic arc r]i meeting ai perpendicularly on both sides 
of ai. Figure [7] shows a lift of ai and r]i in the universal covering. Here ai lifts to ai and 
ai* and r/i to 77^. Note that ai and ai* have the same orientations with respect to r/^. In the 
covering there exist two points, s' G a\ and s* G ai*, on opposite sites of 7/^ and at the same 
distance from r/^, such that s' and s* are mapped to the same point s G 02 by the covering map. 
Observe that we can always find s' and s*, such that the distance r\ from r/^ is equal to . 
We denote by 02' the the geodesic from s' to s* . From we obtain two isometric right-angled 
geodesic triangles. Since intersects s' and s* under the same angle 0, the image 02 of 
under the universal covering map is a smooth simple closed geodesic, which intersects a\ exactly 
once. Hence we can incorporate into our homology basis for S. Applying the cosine formula 
to one of the isometric triangles, we obtain: 



cosh( 



02 ^ 



cosh(ri) cosh( 



^1. 



where r\ 



a\ -ti 



Here the length of rji can be calculated from a decomposition of Qi into a Y-piece (see equation 
(fT4l). We have: 



^^"^^1^ sinh(^) 



cosh(^) _ _^^fVi- 



\ 



'cosh(f) 
sinh(^^) 



+ 1. 



For further calculations we also need the angle 9. From hyperbolic geometry we obtain: 

cos{9) = tanh( "^^^ ) coth(^) 



(12) 
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In Qi there exists likewise a unique shortest geodesic arc rj2 meeting a2 perpendicularly on both 
sides of 02- This arc can be seen in Figure^ Now 02 and ai intersect exactly once under the 
angle 6. Consider the right-angled triangle with sides ^2 and Here r2 contains information 
about the twist parameter t2 with respect to 02- We have: 

Oil 

cos{9) = tanh(r2) coth(— ). 

Together with equation (I12p . we obtain: 

tanh(^)coth(f ) . 2r2 2r2 

tanh(r2) = and t2 = mm( , 1 ). 

coth(i|-) 02 a2 

We will now look for a suitable ai2- Consider again the lifts of oi in Figure^ Consider the two 
points, S' G Qi' and S* £ ai* on the opposite side of s' and s* with respect to the intersection 
point with rji and with distance oi — ri from r][. S' and 5* are mapped to the same point S G Q 
by the covering map. Connecting these points we obtain a geodesic arc a'^2) which maps to a 
scg ai2 in Q. It follows from the intersection properties with ai and 02 that ai2 is in the free 
liomotopy class 01(02)"^- Using the same reasoning as for 02, we can find it's length and the 
twist parameter ti2- 



4.1 Upper bounds based on trirectangles 

Let C be a hyperbolic cylinder with baseline 7. Let T C C be a hyperbolic trirectangle, which 
has one side on 7 C C and such that the two adjacent sides meet d perpendicularly. A lift of 
C together with T into the universal covering is shown in Figure O 




Figure 5: A lift of a cylinder C with baseline 7 with a tube T{j) containing the trirectangle T 
in the universal covering. 

We assume that T is given in Fermi coordinates by 

T = ^l^{{{t,s) \se[0,a2{t)],te [0,d]}), 

where, by abuse of notation, d denotes the length of d. We assume furthermore that the shortest 
side of r, which is perpendicular to d has length w and intersects d in ^(0, 0). For t G [0, d], we 
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cosh(t) \ . , . , / 1 



have 

02(0 = arctanh ( ) \ ] , where d < arcsinh ( — — — - ) . (13) 

\coth(t(;) / \smh(?i;) / 

This follows from the geometry of trirectangles. With Theorem 13. H we obtain the following 
lemma: 

Lemma 4.1 Let C be a hypberbolic cylinder with baseline 7. Let T{'y) C C be a topological tube 
around 7, whose boundary segments consist of either boundary arcs of trirectangles or curves 
contained in a distance set dZrij) for some r. Then we have for H(s) = 2 arctan(exp(s)) and 
Q^it) = '-^\so=a.it)■a[{t): 

'^'^^ - qi{tf+qi(t)q2{t)+q2(tf 



1 + 



dt > cap(r(7)) > j — — ^^^^ " — — dt 



H{a2{t))-H{ai{t)) " ^'^^-y H{a2{t)) - H{a,{t)) 



// a boundary segment dT = '0({(i,-s) | t G [0,d],s = a2(t),0 < w = 02(0) < a2{d)}) C dT{'y) 
belongs to a trirectangle T then 



, / /coth(w) + cosh(t) \ , , , sinh(t) 

H{a2{t)) = 2arctan W and q2{t) - ^ ' 



coth(w) - cosh(t) j y^coth(u;)2 - cosh(t)2 



Such a tube T(7) is depicted in Figure\^ This lemma follows directly from equation (jl3p . the 
formulas for H{s) and qi{t) and the fact that cosh(arctanh(a)) ~ \/l — 2;^ and exp(arctanh(a;)) = 

^J^^■ A parametrisation of other boundary parts of T{'^) can be easily obtained by adapting 
the above formulas. 

The relation with our problem of providing estimates for the energy of dual harmonic forms is 
the following. Let S" be a Riemann surface and A a canonical homology basis. If 7 = aT-(j) for 
some aT-(j) G A and T(aT-(j)) is contained in •S'q^j.j C C, then the capacity of r(a,-(j)) provides 
an upper bound for the entry qa of the period Gram matrix Qs- 

cap(r(a^(j))) > qu. 

We will make use of these formulas in the following section. 

4.2 Upper bounds for the entries of Qs based on the geometry of Q-pieces 

We will establish estimates for all entries of the period Gram matrix based on the geometry 
of the Q-pieces (Qi)i=i,3,..,2g-i- Following the approach given in section 3.1, it is sufficient to 
construct suitable functions on 

S^r\Qi, where 7 G {a^, a^(i), aj^(j)}, for i G {1, 3, .., 2^ - 1}. 

In this section we will only show how to obtain estimates for qu based on the geometry of 
Qi. These estimates will only depend on the Fenchel-Nielsen coordinates (/3i,a2,i2) of Qi. In 
the same way we obtain estimates for (722 and qi2 based on the coordinates, (/3i,ai,ti) and 
ai2; ^12) of Qi- Proceeding the same way on the remaining Q-pieces and combining these 
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estimates as described in section 3.2 we finally obtain estimates for all entries of the period 
matrix. 

To obtain an upper bound for qn we embed H Qi into a hyperbolic cylinder C with baseline 
Q!2 and denote this embedding by the same name. To obtain an estimate on qn we will give a 
parametrization of 

Sa2 n Qi c c 

based on a decomposition into trirectangles. We will give a description of H Qi that depends 
only on the length of 02, /3 = /3i and the twist parameter t2- 

To obtain this parametrization, we first cut open Qi along 02 to obtain the Y-piece Yi with 
boundary geodesies f3, 02' and 02" ■ 02' and 02" have length 02 (see Figure\^. Denote by b the 
shortest geodesic arc connecting 02' and 02" ■ Now we cut open Yi along the shortest geodesic 
arcs connecting /3 and the other two boundary geodesies. We call the surface, which we obtain 
by cutting open Yi along these lines. By abuse of notation, we denote the geodesic arcs in Y( by 
the same letter as in Yi. The geodesic arc b divides Y^ into two isometric hexagons Hi and H2- 
This decomposition is also shown in FigurelSi 




Figure 6: Decomposition of Yi into isometric hexagons Hi and H2- 

In Hi b is the boundary geodesic connecting ^ and Denote by 6' the shortest geodesic arc 
in Hi connecting b and the side opposite of b of length ^. By abuse of notation we denote this 
side by ^. We denote by 6" the arc in H2 corresponding to 6' in Hi. Let 6 be the geodesic arc 
in Y(, such that 6 = 5' U 6" . It is easy to see that the geodesic arc in Qi that maps to 6 in Y( 
constitutes the intersection of the cut locus of 02 with Qi. We denote this geodesic arc in Qi 
and Yi also by 5: 

6 = CL{a2)nQi. 

We denote furthermore by a' the geodesic arcs connecting ^ and ^ in Hi and by a" the corre- 
sponding arc in H2- Let a be the length of such an arc. 

6' divides Hi into two isometric right-angled pentagons Pi and P2. Let Pi be the pentagon that 
has ^ as a boundary. To establish the parametrization for H Qi, we divide Pi into two 



17 



trirectangles. 

Let c be the geodesic arc in Pi that emanates from the vertex, where ^ and 6 intersect and that 

meets ^ perpendicularly. It divides ^ into two parts, a' and a" (see Figure\^. c divides Pi 

into two trirectangles, Ti, which has the bomidary a' and T2, which has the boundary a" . 

To obtain an upper bound for qn, we have to know the geometry of Ti and T2. To apply the 

results from the previous subsection, we have to know the lengths a, a', a" and |. 

In the following subsection we will also need the length of c, which we will calculate here. To 

obtain these lengths, we will use the geometry of Hi, Pi, Ti and T2. 



From the geometry of the hyperbolic pentagon Pi we have: 

1.^^^ • v,^"2^ . .fb . b cosh(|) 

cosh(-) = smh(y)smh(-) 4^ smh(-) = — (14) 

cosh(5') = sinh(— ) sinh(a). (15) 

Hence we obtain b from equation ([HI). 

From the geometry of the hyperbolic hexagon Hi and as sinh(6) = 2 sinh(|) cosh(|) and cosh(6) = 
2cosh(|)2 — 1, we obtain now a: 

sinh(fe) cosh(^) b 02. 

coth a = — — — = tanh - cosh — . 16) 

1 + cosh(oj 2 2 

From equation ([TS]) and (fT6|) and as sinh(arccoth(x)) = ^J^ ^ we obtain 5' = |: 

cosh(y) = ^ ""^(^^ (17) 
/tanh(|)2cosh(^^)2 - 1 

Finally from the geometry of hyperbolic trirectangles Ti and T2 we obtain : 

, ^ „^ cosh(l) ^ ^ 

coth(a" = riTK- 18 

^ ' tanh((5') ^ ^ 

sinh(c) = cosh((5') sinh(^). (19) 

As a' = ^ — a" , it remains to express a" in terms of a2 and b. From equation (|17p and (jlSp we 
obtain with 7 — w — - — r-r-rr = / that 

tanh(arccosn(a;jJ y/x^ — l 



coth(Q;") 



cosh(|)2sinh(^^) 
cosh(|)2 sinh(^a)2 _ (sinh(|)2 cosh(^a)2 _ cosh(|)2) 



"J^^^^n 2 ^"^^"^^^2 

\2 _ ^^„t,,'^\2 



As sinh(x) = cosh(x) — 1, this simplifies to 



■ b.2, 1 /«2^ 



coth(a") = cosh(-)nanh(y). (20) 
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Finally we obtain for c from equation (|14p . (|17p and (I19|) and as sinli(x)^ = cosh(3;)^ — 1 



sinh(c) = (21) 
tanh(|)2cosh(^)2 - 1 

Hence can express the length of 6, a, q" and a' = ^ — a" in terms of and /3 - see equation 
(fT4l) . (fTBI) and ([20]), respectively. 



With these formulas we are now able to obtain a description of the boundary of S^^ H Qi C C. 
Consider now b C Y[. 6 divides into two isometric hexagons. Let H[ be the hexagon that 
contains as a boundary geodesic and H2 be the hexagon that contains 02" as a boundary 
geodesic 5 forms the cut locus of 02 in Qi. Denote by Q\ the surface that we obtain if we cut 
open Qi along S. Q'l is a topological cylinder around 02- A lift of Q'l in the universal covering 
is depicted in Figure Q 





Figure 7: Lift of Q'^ into the universal covering. 

By abuse of notation we denote by H[ and H2 two hexagons in this lift, which are isometric to 
the hexagons with the same name in and that are adjacent along the lift d2 of 02- We keep 
the notation from Y(, but denote by 5i C H- , for i G {1,2}, the two sides corresponding to 5 in 
Y^. In the lift of Q\ the two hexagons H[ and H2 are shifted against each other by the length 
t2 ■ OL2- It is easy to see from Figure^ how to parametrize Sa^ H Qi in a cylinder C around 02- 
Here all boundaries are boundaries of trirectangles, which are isometric to either T\ or T2. Hence 
n Qi = T(a.2) satisfies the conditions of Lemma 14.11 Adapting the formulas given in this 
lemma, we can find an upper bound for 

cap(S'o2 n Qi) > qx\. 

Remark: From numerical simulations we get the following impression. If the length of ol or a" 
is near the maximal possible value arcsinh( ^.^^^^ ) or arcsinh(^r^^j^), respectively, then our test 

function is not well suited for the capacity problem and overestimates the capacity of S^^ H Qi 
and therefore also E{a\). Here we can improve our estimate by cutting of these pointy edges at 
the endpoints of b\ and by replacing them with a suitable boundary of constant distance from 
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02 (see Figurel^. 

We also obtain an upper bound for E{ai) from the capacity cap(C(a2)) of a collar around 02 (see 
equation ([2])). The upper bound from cap(S'Q2 Pi Qi) seems to be improving the estimate from 
cap(C(a2)) if min{a, |} is small. For min{a, |} > 4 hardly any improvement can be obtained. 

4.3 Lower bounds on the entries of Qs based on Q-pieces 

Consider the primitive Fi of ai in Q\ = D Qi C C . The two geodesic arcs 5i and 62 
corresponding to 6 C Qi constitute C L{a2y^'^ r\ dQ'i- We will use the theoretical approach from 
section 3 to obtain a concrete lower bound for 

qu = Es{Fi) > Eb{Fi), where B = B2 = Sl""^ H Q[. 

We will give a suitable construction for B = 82^'^ n Q'^ in the following. To this end we lift 
Q']^ into the universal covering as in the previous subsection (see Figure Q • We use the same 
notation for the geodesic arcs that occur. The important cutout from Figure is depicted in 
Figure 




Figure 8: The area B and the construction of skewed Fermi coordinates . 

Let B be the hatched subset in the lift of Q'^ in Figure [21 We will give an exact description and 
parametrization of B in the following. 

The boundary of B contains the lines 5i and 82- For each point pi £ 61 there exists a point 
P2 S such that pi and p2 map to the same point p on 6 <Z Qi- We may assume WLOG that 

Fi{p2) - Fi{pi) = 1 for all pi G 61 

We will describe S as a set of lines, where each line Ip connects pi and p2- Ip can be described 
in the following way. From pi we go along the geodesic that meets 02 perpendicularly until we 
meet dZb{d2)- We call this intersection point p'^ and the geodesic arc that forms 7^. Let P2 be 

the point on dZb{a2) on the other side of 02 that can be reached analogously, starting from p2. 

2 

We now go along the geodesic arc that connects pi and P2. We call this arc 7^. Then from P2, 
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we move along the geodesic arc connecting ^-nd p2- We call this arc 7^. We define Ip as the 
line traversed in this way. Let B be the disjoint union of these lines: 

B=\S{lp} 

pes 

We will use a bijective parametrization of B. We want ay?: (t, s) 1— )• (p{t,s) such that for a 
fixed to £ 92(^0; ■) parametrizes the line Ip that traverses 02 in a point with directed 

distance to from m. We parametrize the sets |J {7p} and |J {7^} in Fermi coordinates with 

p£S pes 

baseline d2- The proper parametrization can be deduced from the geometry of the trirectangle 
T2. 



We will parametrize Zb{a2) D B = IJ {TpI using skewed Fermi coordinates with angle u 

2 pes 

and baseline d2- These are defined in the same way as the usual Fermi coordinates tp (see section 
1), but instead of moving along geodesies emanating perpendicularly from the baseline, we move 
along geodesies that meet the baseline under the angle v. We will not give these coordinates 
explicitly, but will derive the essential information from the usual Fermi coordinates ij). 
Therefore let A be the geodesic arc connecting the midpoints of 5i and 62- The midpoint m of 
A and the endpoints of | are the vertices of a right-angled triangle D (see Figure\^. In our 
case the angle v for the coordinates tp'^ is the angle of D at the midpoint m. It follows from the 
geometry of right-angled triangles that 

cosh(^) = cosh(^)cosh(^), (22) 

where we assume WLOG that the twist parameter t2 is in the interval [0, \]. Otherwise the 
situation is symmetric to the depicted one. We have: 

sinh(- ' — 



2 sin(i/) 

The equations for D imply with sinh(x)^ = cosh(2;)^ — 1: 

sinh(l) ^ ^ 

sin(i/) = (23) 

^cosh(|)2cosh(^^^)2 - 1 

Consider the following geodesic arcs in Zb{d2) H B. For a n G N let a'„ be a geodesic arc of 

length on 0.2 with midpoint m. A intersects in m under the angle v. This is depicted in 
Figure [3 

Let Vj' be a geodesic intersecting A perpendicularly in m. Let furthermore //i and /U2 be two 

geodesic arcs with endpoints on Zb{d2) that intersect rj' perpendicularly, such that each of the 

2 

arcs passes through an endpoint of on each side of A. Let furthermore rj be the geodesic arc 
on rj' with endpoints on /ii and ^12- For fixed n G N, we denote by the length of rj and by jjJ^ 
the length of jii and ji2: 

Tjn = Kv) and = l{jii) = l{ji2)- 
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By choosing usual Fermi coordinates with basehne ry, we can parametrize the strip, whose bound- 
ary hues are ^ui and ^2 and two segments of dZb{a2) (see Figure\^. 

2 

n such strips can be ahgned next to each other to obtain a parametrization of Zb{a2) H B. For 

2 

n — )• 00 we obtain a parametrization ij)'^ of Zb{di2) H B. We have: 

2 

hm n ■ rjn = sin(i/)2a" and hm /x" = A. 

n— >oo n— >oo 

Combining the parametrizations for the several pieces oi B, we may assume that we have a 
parametrization ip that satisfies our conditions. For practical purposes, we extend the parametriza- 
tion if to the geodesies meeting dZb n B perpendicularly in the direction opposite of 0.2- 

2 

Consider a point pi = (p{to,—x) £ 61 and p2 = (p{to,x) £ 82 ■ The function Fi satisfies the 
boundary conditions Fi{p2) = 1 + c and Fi{p2) = c, where c is a constant. As we will see in the 
following, the constant c is not important for our estimate and we assume that c = 0. 
We consider the strip V, where 

V = (p{[to — e,to] X [— where e > if to < and e < if to > 0. 

We will show, how to obtain a lower bound for the energy of -Fi|\/) Ev{Fi) for sufficiently small 
e. We can align these strips to obtain a lower bound for Eb{Fi) < E{Fi). We derive a lower 
bound for the energy of -Fi|v, assuming that 

^ilv>([to-e,io]x{-x}) = Fi{pi)={) and = i^ife) = 1 and 

^il^([io-Mo]x{^}) = ^i(K)=«'i and Fi|^([,^_^^,^]^^.^)=Fi(p'2)=a'2. 



Consider the subset V' of V given by 

V' = ^{[to-eM X [^,^])- 
Let = /tf, otjj^ he a, function defined on V that realizes the minimum 

min{ j \\p^{Df)\\l I / e Lip(y'),/l^([to_,,to]x{^} = «i and /|^([,^_, ^^j^i = 03} (see eq. ®). 
v 

By considering skewed Fermi coordinates as a limit case of Fermi coordinates with respect to an 
imaginary baseline r] and by applying the calculus of variations (see [Ge|, p. 14-16) to the last 
integral in equation ([1]) we obtain that ft^ is given by 



02 -ai , , aiH{^)-a2H{-^) 

If nit, S) = T 7— His) H T r— ^ 

The energy Ev'{p,p{DFtp)) is 



V (Pi/3( io)) 2(arctan(exp(^)) — arctan(exp(— ■!))) 1) I I ( ) 

We can extend Ft^ to a function on V that satisfies the boundary conditions 

Fto\ip{[to-eM^{±x)) = -^l|(^([to-e,to]x{±a:})- 
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Again we choose Ff^, such that it minimizes Ey\^yi {p^{D{-))) with the given boundary conditions. 
We have with Ev\v'{Pip{DFto)) = Ey\y,{Fta): 

Ey\y.{F^,) = ^. ^ f^,^l~""^y, = ^2(to)(a? + (1 - a^f^, (25) 

2(arctan(exp(xjj — arctan(exp(2 jjj 

where x = x(to) depends on t^. For oi = a']^ = Fi{p't^ and 02 = a2 = -Fi(p2), we have by 
construction Ey^Fi) > Ey {pyy{D Ff^^)) . 

Though we do not know the values a'^ and ag, we obtain a lower bound of the energy of Fi on 
V, if we determine the values FtQ^p'i) = ci = ci(to) and Ftg{p2) = C2 = C2(to)) respectively, such 
that these values are minimizing the total energy Ey^p^iFt^^)). As the two arcs 7^ and 7^ have 
the same length, we have to solve the following problem: 
Find ci, C2, such that 1 — C2 = ci <^=> (c2 — ci) = 1 — 2ci, such that 

Ey{p^{DFt,)) = Ey>{p^{DFt,))+Ey\y,{Ft,) 

is minimal. We obtain from equation (|24p and (|25p that ci = ^3(^0)^-2^1 " cover i? with 

a set of such strips Vjq = V , such that these intersect only on the boundary and combine the 
FtglvtQ to a function fi on B. We have Ey{Fi) > Ey{p^{Dfi)). 

As we consider only the energy EB{P(fi{Dfi)) of the projection, the approximation is true in 
the limit case, where e — )• 0. Due to the symmetry of the area B we obtain in total with ci in 
equation ([MI) and ([25]) : 

a" 

m = £(F,) > = 2 / jj|2(|_<i, (26) 

i?y\y/(p<^(Di<t(j)) is monotonously decreasing if the length is x — | increasing. Hence we find a 
simpler approximation for E{Fi), setting x = c (see (|2T]) ) in equation (I25p . In this case, we can 
set k2{t) = k2{ct") for all t. We obtain furthermore a simplified upper bound, if we define our 
test function only on Z^/(a2), where w' = min{a, |}. This way we obtain: 

>E[Fi)-qii>-————. (27) 



2(arctan(e'"') -arctan(e""'')) ~ ^ ' k2{a") + 2ki 

The lower bound that depends only on a2,t2 and /3i can be obtained by expressing | and A and 
1/ in terms of these variables (see equation (jl4p , (j22p and (|23p ) . Using the parametrization of T2 
in equation (I25p we can express the lower bound in terms of the variables 02 , ^2 and /3i . This 
way we obtain explicit values in equation (j26p . 



Remark: From numerical simulations we get the impression that the lower bound for E{Fi) 
is only usable, if the twist parameter t2 is small. The reason for this seems to be that using 
the projection strongly underestimates the energy. The result gets worse with increasing t2 and 
increasing length of a2- 

4.4 Summary and discussion 

Let Qi be a Q-pieces from our decompositon. We use the notation as in the previous subsections. 
For j G {i,T{i),iT{i)} let Bj be the set corresponding to B = B2 = -6(02) in the previous 
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subsection. Let fT(j) defined on Bj analogously to the function fi = /r(2) on B2 (see equation 
(j26p . Let furthermore /ir(j) on -BjT-(i) be the function corresponding to fi on B2 which we use 
to obtain a lower bound for E{ai + on Sa.^^.^ . We obtain the following upper and lower 

bounds for the entries of the matrix Qs, where Si = Sa^: 

cap{Sr(i) n Qi) > qu > EB^(i){p^{D fi)). 

Note that by definition Qi = Qr(i) is the Q-piece containing Ui and aT-(j). Furthermore 

gi,r(i) < ^ (cap(5i^^(i) n Qi) - EB^MDf^i^i)) - Eb^^,^ {p^{Dfi))^ and 
qi,r{i) > \ (-E^B^.^wlPvl-^/i.rCi))) - cap(S'i n Qi) - cap(5^(j) n Qi)) . 
For the entries gj / involving estimates from different Q-pieces Qi and Qi we have (see equation 

cnD: 

< ^ (cap(S'^(i) n Qi) + cap(5^(i) n Qi) - Eb^^^-^{p^{D fi) - EB^f^,){p^{D fi))^ . 

From numerical tests we get the impression that our estimates for an entry qr{i)^T{i) good, if a 
large part of the cut locus of CL{ai) is contained in Qi and if Oiti is small. The first condition 
is fulfilled, if the boundary geodesic /3j is small. The estimate for qu involving different Q-pieces 
can only be sharp, if qu is small. Again, if /3j and /3; are small and if both and are small, 
we obtain a good estimate. The methods seem to work best, if all ctj,j € {i, r(i), ir(i)} are 
small. This justifies the choice of ai in equation (|lip . Note that by |BSel| it exists a canonical 
homology basis for a Riemann surface of genus g, where the largest element is of order g. Hence 
this condition can in prinpicle be satisfied for small g. 

The advantage of the method is that information about the geometry of the surface can be 
incorporated. Suppose, for example, that the geometry of a Y-piece Yi attached to the Q- 
piece Qi from our decompositon is known. Then the the cut locus CL{a2) H {Qi U Yi) can be 
calculated. Incorporating this information we obtain better estimates for E{ai). Information 
about isometries of the surface can also be incorporated. This was used in Example 2. 
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